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ABSTRACT 

The Sunyaev Zel'dovich effect (SZ effect) is a complete probe of ionized baryons, the 
majority of which are likely hiding in the intergalactic medium. We ran a 512 3 ACDM 
simulation using a moving mesh hydro code to compute the statistics of the thermal 
and kinetic SZ effect such as the power spectra and measures of non-Gaussianity. The 
thermal SZ power spectrum has a very broad peak at multipole I ~ 2000 — 10 4 with 
temperature fluctuations AT ~ 15/^K. The power spectrum is consistent with available 
observations and suggests a high o"s ~ 1.0 and a possible role of non-gravitational 
heating. The non-Gaussianity is significant and increases the cosmic variance of the 
power spectrum by a factor of ~ 5 for I < 6000. 

We explore optimal driftscan survey strategies for the AMIBA CMB interferometer 
and their dependence on cosmology. For SZ power spectrum estimation, we find that 
the optimal sky coverage for a 1000 hours of integration time is several hundred square 
degrees. One achieves an accuracy better than 40% in the SZ measurement of power 
spectrum and an accuracy better than 20% in the cross correlation with Sloan galaxies 
for 2000 < I < 5000. For cluster searches, the optimal scan rate is around 280 hours 
per square degree with a cluster detection rate 1 every 7 hours, allowing for a false 
positive rate of 20% and better than 30% accuracy in the cluster SZ distribution function 
measurement. 

Subject headings: Cosmology-theory-simulation-observation: SZ effect, cosmic microwave 
background, large scale structure, cluster. 
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1. Introduction 

Big bang nucleosynthesis (BBN) and cosmic microwave background (CMB) experiments such 
as Boomerang (Netterfiled et al. 2001) and DASI (Pryke et al. 2001) predict that the ordinary 
baryonic matter accounts for about 5% of the total matter in the universe. But the directly 
observed components, such as stars, interstellar medium and intracluster gas, only account for a 
few percent of this baryon budget, while more than 90% baryons have escaped direct detections 
(Persic &; Salucci 1992; Fukugita, Hogan & Peebles 1998). This is known as the missing baryon 
problem. The missing baryons are believed to be in the form of the intergalactic medium (IGM) 
and have been difficult to detect directly. To understand their state, such as density, temperature, 
peculiar velocity and metalicity, stands as a major challenge to both observation and theory, and 
is crucial to understand the thermal history of the universe and galaxy formation. 

The IGM has various direct observable tracers. (1) Neutral hydrogen absorbs background 
light of quasars and produces the Lyman-a forest. (2) Ionized electrons have thermal and peculiar 
motions and are capable of scattering CMB photons and generate secondary CMB temperature 
fluctuations, which are known as the thermal and kinetic Sunyaev Zel'dovich effects (SZ effect), 
respectively. Their precision measurements are becoming routinely available with the devoted 
CMB experiments such as AMIBA (Array for Microwave Background Anisotropy 2003) and SZA 
(Sunyaev Zel'dovich array). (3) Ionized electrons and protons interact with each other and emit 
X-rays through thermal bremsstrahlung and contribute to the soft X-ray (0.5 — 2 keV) background 
(XRB). Several other tracers have been proposed, including X-ray absorption techniques (Perna & 
Loeb 1998), but they depend strongly on chemical compositions, and would be difficult to associate 
with direct IGM properties. These tracers probe different IGM phases and help to extract the IGM 
state. For example, the XRB flux upper limits place constrain on the gas clumpyness and suggest 
a potentially strong role of feedback (Pen 1999). 

Among these tracers, the SZ effect is a particularly powerful IGM probe. (1) It provides a 
complete sample of the intergalactic medium. All free electrons participate in Thomson scattering 
and contribute to the SZ effect. The Thomson optical depth from the epoch of reionization z ~ 10 
to the present is r ~ 0.1, which means that about 10% of CMB photons have been scattered by 
electrons. Since the number of CMB photons is much larger than the baryon number and the 
ionization fraction of our universe is larger than 99%, most baryons contribute to the SZ effect. 
Compton scattering does not depend on redshift and is not affected by distance or the expansion 
of the universe. So, the SZ effect can probe the distant universe. Furthermore, the SZ effect 
does not strongly depend on gas density and probes a large dynamic range of baryon fluctuations. 
In contrast, the X-ray emission measure depends on density squared, and primarily probes the 
densest IGM regions at low redshift. The Lyman-a forest probes the neutral IGM, which in 
ionization equilibrium depends also on the square of the baryon density. The neutral fractions only 
accounts for a tiny fraction of total baryons. Extrapolations over many orders of magnitude are 
required to understand the bulk of the baryons. (2) It is straightforward to understand. The simple 
dependence of the SZ effect on the ionized gas pressure or momentum does not put as stringent 
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requirement on simulation resolution, nor as accurate an understanding of the gas state such as 
metalicity, temperature, velocity and ionization equilibrium, as X-rays and Lyman-a properties 
do. Pressure is the total thermal energy content of the gas, and finite volume flux conservative 
simulations are particular amenable to its modelling. It is also easier to model analytically. (3) It 
has strong observational potential. In this paper we will primarily consider the thermal SZ effect, 
which is easier to observe. Its unique dependence on frequency allows us in principle to disentangle 
the SZ effect from the contamination of the primary CMB and other noise sources (Cooray, Hu k 
Tegmark 2000). Hereafter, if not otherwise specified, the SZ effect always means the thermal SZ 
effect. 

In SZ observations, all redshifts are entangled and smear some key information of the inter- 
vening IGM: its spatial distribution and time evolution. These properties around redshift z ~ 1 are 
sensitive to many cosmological parameters and thus potentially promising to break degeneracies 
in cosmologies from the primary CMB experiments. One can hope to resolve the equation of the 
state of the dark energy. We are currently performing simulations for different cosmologies, which 
are degenerate in the primary CMB, to test the potential of such a procedure (Zhang et al. 2002). 
If combined with other observations, the SZ effect will become even more powerful. For example, 
with the aid of a galaxy photometric redshift survey, the evolution of the 3D gas pressure power 
spectrum and its cross correlation with the galaxy distribution can be extracted (Zhang k Pen 

2001) . 

The above analysis depends on a detailed quantitative theoretical understanding of the SZ 
statistics. Analytical approaches considered in the past strongly depend on various assumptions. 
The Press-Schechter approach as adopted by Cole k Kaiser (1988); Makino k Suto (1993); Atrio- 
Barandela k Mucket (1999); Komatsu k Kitayama (1999); Cooray (2000); Molnar k Birkinshaw 
(2000) requires ad hoc models for the gas profile in halos, whose shape and evolution are uncertain. 
The hierarchical method as proposed by Zhang k Pen (2001) strongly depends on the gas-dark 
matter correlation, which is also poorly understood. These estimates can be significantly improved 
by high-resolution simulations. But past simulation results (Scaramella et al. 1993; da Silva et al. 
2000; Refregier et al. 2000; Seljak, Burwell k Pen 2001; Springel et al. 2001) disagreed on both 
amplitude and shape of the SZ power spectrum (see Springel et al. (2001) for a recent discussion). 
Simulation resolution may play a key role in this discrepancy since the SZ effect is sensitive to small 
structures, as discussed by Seljak, Burwell k Pen (2001). Differences between code algorithms may 
also cause part of the discrepancy. To address these problems, we ran the largest adaptive SZ 
simulation to date, a direct 512 3 moving mesh hydro (MMH) code (Pen 1998a) simulation, to 
reinvestigate the SZ statistics. For a better understanding of numerical and code issues in the SZ 
effect, our group is carrying out a series of simulations with identical cosmological parameters and 
initial condition but using different codes with different resolutions (Codes Comparison Program 

2002) . 

Apart from simulation issues, a problem that has received little attention is analysis strategies 
for SZ data. Several interferometric arrays such as AMIBA and SZA are under construction, but 
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detailed data analysis models are still in their infancy. With our simulated SZ maps, we can 
estimate the sensitivity and accuracy of SZ observations. Given a sky scan rate, how accurately 
can the SZ power spectrum be measured? How many cluster can be found? How accurately can 
the SZ decrement be determined for individual clusters? What are optimal strategies for these 
measurements? In this paper, we will take AMIBA as our target to address these questions. 

This paper is organized as follows: in §2, we describe the SZ effect, its statistics and our method 
to analyze these statistics. We then present our simulation results of these statistics and possible 
constrain from observations (§3). In §4 we simulate AMIBA driftscan observations to estimate the 
accuracy of AMIBA measurement of these statistics. We conclude in §5. 



2. Formulation and Background Review 

The thermal SZ effect causes CMB temperature and intensity fluctuations in the sky. At a 
given position given by a unit vector n pointing from the earth to some point on the sky, the 
fractional change in the CMB intensity SI U /I U depends on both the observing frequency v and the 
integral of the gas pressure along the line of sight (Zel'dovich and Sunyaev 1969): 

^ = -2y(n)Sj[x(u)]. (1) 

For the CMB intensity, the spectral dependence for inverse compton scattering off non-relativistic 
electrons is described by Sj(x) = -0^1 (2 — x/[2 tanh(x/2)]) where x = hv / (/cbIcmb) = ^/57 Ghz. 
The corresponding fractional change in CMB temperature 6(n) = A7 ^ v ^ TO - > is then: 

9(n, v) = -2y{n)S T [x(v)] = -2y (2 - x/[2 tanh(x/2)]) . (2) 

In the Rayleigh-Jeans limit (x <C 1), Si(x) = St(x) = 1. Si(x) reaches the highest response, 
S(x) ~ 1.6 at v ~ 100 Ghz, which is well matched to the AMIBA frequency range of 80-100 Ghz. 
St(x) monotonously decreases with frequency. We show Si(x) and St(x) in fig. (1) where the 
frequency responses of various experiments which already have the CMB observation data at small 
angular scales (Z > 2000) such as ATCA (Subrahmanyan et al. 2000), BIMA (Dawson et al. 2001), 
CBI 1 , SUZIE (Church et al. 1997) and VLA (Partridge et al. 1997) are shown. The gas pressure 
dependence is described by the Compton y parameter 

a T [ lre cf t [ lie 

y(n) = o / n e (lh)kBT e (lh)dl = » / PJlnjdl. (3) 

m e c 2 J ■ m e c 2 J 

T e , n e and P e are the temperature, number density and pressure of free electrons, respectively, dl 
is the proper distance interval along the path of CMB photons, or, m e and c have their usual 
meanings as the Thomson cross section, electron mass and the speed of light. Since only free 
electrons participate in Thomson scattering of the CMB photons, the integral is cut off at the 
epoch of reionization Z rc . 



1 http://www. astro. caltech.edu/ tjp/CBI/ 
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Fig. 1. — The SZ spectral dependence of CMB brightness fluctuation Si(v) and temperature 
fluctuation St{v)- 5I U /I^ MB = —2ySi(u). AT/Tqmb = — 2ySr(i/). We show the responses of 
various experiments which are capable of detecting SZ signals such as AMIBA: 90 Ghz (Array for 
Microwave Background Anisotropy 2003), ATCA: 8.7 Ghz (Subrahmanyan et al. 2000), BIMA: 
28.5 Ghz (Dawson et al. 2001), CBI: 26-36 Ghz, SUZIE: 142 Ghz, 217 Ghz, etc. (Church et al. 
1997), SZA: 26-36 Ghz and 85-115 Ghz and VLA: 8.4 Ghz (Partridge et al. 1997). S(u) -» 1 when 
v <C 50 Ghz. At 217 Ghz, there is no thermal SZ effect, enabling the measurement of other CMB 
anisotropies such as the kinetic SZ effect. 
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Some important statistics of the SZ effect immediately come to mind. The first order quantity 
is y, the mean y parameter averaged over the whole sky, which measures the total thermal energy 
content of the universe. The angular variation in y can be parameterized by the two point correlation 
function of the temperature fluctuation 6, or equivalently the angular power spectrum C\. For a 
Gaussian random field, these two parameters would describe the statistics completely. Since the SZ 
effect is dominated by non-linear structures, non-Gaussianity may be significant. So we investigate 
higher order statistics such as the skewness and kurtosis of the y parameter to quantify the non- 
Gaussianity. 

The SZ effect contains contributions from all redshifts and it is challenging to recover the 
smeared redshift information. We have shown that cross correlating the SZ effect with a galaxy 
photometric redshift survey, we can infer the redshift resolved IGM pressure-galaxy cross correlation 
and the IGM pressure auto correlation (Zhang & Pen 2001). This method is robust, but does not 
capture all the information in the SZ observation. In this paper, we utilize the one point distribution 
function (PDF) of the y parameter and the distribution of peaks in y to extract more information. 

Since smoothing is always present for a real experiment with a finite beam, we calculate the 
statistics of the y parameter smoothed on a given angular scale 8. If we are interested in virialized 
objects, for example clusters and groups of galaxies, we would expect them to be peaks in the 
smoothed or filtered y maps. N{y > y p ), the cumulative distribution function (CDF) of peaks with 
smoothed y parameter bigger than certain value y p , is the raw observable. The y p CDF is the SZ 
analog to a luminosity function. If we choose a top hat window so that the observation cone is large 
enough to include an entire object such as a cluster and is small enough that typically no more than 
one such object can be found along each cone, then when the cone is centered at the center of each 
object, a peak y = y p appears in the smoothed map. This y p is directly related to the total gas mass 
and temperature of individual object. Assuming a halo to be isothermal, the total mass M of a halo 
is related to the gas temperature T by M/M 8 = (T/T 8 ) 3 / 2 . M 8 = 1.8 x 1O 14 (O o /O.3)/i" 1 M is the 
mass contained in a 8h _1 Mpc sphere of the universe of mean density today, which is roughly the 
mass scale of clusters. T$(z) is the corresponding temperature of a halo with mass Mg at redshift 
z. We then obtain 

= Msf l ^^knTs(M_\ 5/3 () 
yp nm H d 2 A Afl m e c 2 \M 8 J ' K) 
Here, Af2 is the solid angle of the cone, f g is the gas fraction of halos and cIa is the angular diameter 
distance. For clusters and groups, the typical angular size at z = 1 is about l'. For the present 
cluster number density n(T > 2keV) ~ 10~ 5 /i -3 Mpc 3 (Pen 1998b), the average number of clusters 
in a cone with angular radius 8 ~ 20' projected to z ~ 2 is about one allowing for the evolution of 
cluster number density. So, the size of the smoothing scale for the peak analysis should be between 
these two scales. In this case, the N(y > y p ) is just the number of halos with y > y p , which is given 
by: 

fro fOO i 

N(y >y P )= X 2 dr/Vl - Kr2 / — dM. (5) 

JO JM{y p ,z) & M 

■j^j(M,z) is the halo comoving number density distribution function and is well described by the 
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Press-Schechter formalism (Press & Schechter 1974). x, r and K are the comoving distance, radial 
coordinate and curvature of the universe, respectively. The subscript 're' means the reionization 
epoch. M(y p ,z) is the mass of the halo with smoothed y = y p given by eqn. (4). Eqn. (5) has 
two applications. Firstly, given a SZ survey and the best constrained cosmological parameters 
as determined by CMB experiments, Type I a SN, weak lensing, etc. and a SZ survey, the only 
unknown variable in eqn. (5) is Tg(z)f g (z), which then is uniquely determined. Tg is robustly 
predicted since it is mainly determined by Mg through hydrostatic equilibrium and has only weak 
dependence on the thermal history. For example, comparing the cluster temperature function 
as inferred from simulations with the Press-Schechter formalism, Pen (1998b) found that Tg = 
4.9(1 + z)n 2 /3 n(z)°- 283 keV for a ACDM universe. The above relation is sufficient to extract the 
evolution of the gas fraction f g , which is very sensitive to the thermal history. For example, non- 
gravitational energy injection decreases f g . Secondly, the number of clusters strongly depends 
on cosmology, as characterized by -j^. Given a good understanding of cluster temperature and 
gas fraction, eqn. (5) allows one to constrain cosmological parameters. In combination with the 
measurement of cluster redshifts, this method is more sensitive (Weller, Battye & Kneissl 2001). In 
§4, we will study the statistics of these peaks for maps filtered in optimal ways to measure clusters 
of galaxies. 

With the measurement of the cluster SZ temperature distortion and follow up X-ray ob- 
servations of cluster X-ray flux Fx and X-ray temperature T, cosmological parameters can be 
constrained. For a cluster with electron number density n e , proper size L and angular size 9, 
Fx oc n^K(T)L? / d\ and y oc n e TL. Here, A(T) is the X-ray emissivity temperature dependence. 
Then the luminosity distance (Il(z) oc j^-A^ft a+zp ■ The only uncertainties in this relation are the 
intracluster gas profile and metalicity which affects the X-ray emissivity. These properties could be 
modelled and are potentially observable. In this sense, the SZ effect can be used as a cosmological 
distance indicator (Silk & White 1978; Barbosa et al. 1996; Mason, Myers & Readhead 2001; Fox 
& Pen 2001). Since dark energy dominates at low redshift where the SZ observation and X-ray 
observation of clusters are relatively straightforward, SZ clusters are a potential probe to constrain 
the equation of state for the dark energy. 



3. Simulations 

In this section, we describe our SZ simulation used to investigate the above statistics. We will 
also use our simulation to provide SZ maps for our estimation of the sensitivity of upcoming SZ 
experiments (§4) and its effect on data analysis strategies. We used a moving mesh hydrodynamics 
(MMH) code (Pen 1998a). It features a full curvilinear total- variation-diminishing (TVD) hydro 
code with a curvilinear particle mesh (PM) N-body code on a moving coordinate system. The full 
Euler equations are solved in an explicit flux-conservative form using a second order TVD scheme. 
The curvilinear coordinates used in the code are derived from a gradient of the Cartesian coordinate 
system. If x 1 are the Cartesian coordinates, the curvilinear coordinates are £ l = x l + dl(p(^). The 
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transformation is completely specified by the single potential field t). The potential deformation 
maintains a very regular grid structure in high density regions. The gravity and grid deformation 
equations are solved using a hierarchical multigrid algorithm for linear elliptic equations. These 
are solved in linear time, and are asymptotically faster than the FFT gravity solver. At the same 
time, adaptive dynamic resolution is achieved. During the evolution any one constraint can be 
satisfied by the grid. In our case, we follow the mass field such that the mass per unit grid cell 
remains approximately constant. This gives all the dynamic range advantages of smooth particle 
hydro (SPH) combined with the speed and high resolution of grid algorithms. The explicit time 
integration limits the time step by the Courant condition. To achieve a reasonable run time, we limit 
the compression factor to a factor of 5 in length, corresponding to a factor of 125 in density. Most 
SZ contributions arise below such densities, giving a diminishing return to go to higher compression 
factors. 

The parameters we adopted in our 512 3 simulation are f^o = 0.37, 17 a = 0.63, 17b = 0.05, 
h = 0.7, as = 1.0, power spectrum index n = 1, box size L = 100/i -1 Mpc and smallest grid spacing 
40/i _1 kpc. The simulation used 30 GB memory and took about three weeks (~ 1500 steps) on 
a 32 processor shared memory Alpha GS320 at CITA using Open MP parallelization directives. 
During the simulation we store 2D projections through the 3D box at every light crossing time 
through the box. The projections are made alternatively in the x, y, z directions to minimize the 
repetition of the same structures in the projection. We store projections of thermal SZ, kinetic 
SZ, gas and dark matter densities. For the thermal SZ, we store 2 Ay = ^^P e L as given by 
eqn. (3). Our 2D maps are stored on 2048 2 grids. As tested by Seljak, Burwell & Pen (2001), 
this preserves all the information down to the finest grid spacing. After the simulation, we stack 
the SZ sectional maps separated by the width of simulation box, randomly choosing the center of 
each section and randomly rotating and flipping each section. The periodic boundary condition 
guarantees that there are no discontinuities in any of the maps. We then add these sections onto 
a map of constant angular size. Using different random seeds for the alignments and rotations, we 
make 40 maps of width 1.19 degrees to calculate the SZ statistics. The mean y parameter in these 
40 maps y ~ 4 x 10" 6 . The mean y parameter is still below the upper limit 1.5 x 10 5 from COBE 
FIRAS (Fixsen et al. 1996). One typical thermal SZ sky map and a kinetic SZ sky map at the 
same field of view are shown in fig. 2 and fig. 3, respectively. With these maps, we calculate the 
SZ power spectrum, the SZ non-Gaussianity, the y PDF and y v CDF. 



3.1. The SZ power spectrum 

The power spectrum in the Rayleigh-Jeans regime averaged over 40 maps is shown in fig. 4. 
The thermal SZ power spectrum is fairly flat, with a broad peak at / ~ 2000-10 4 and a typical 
fluctuation amplitude O ~ 5.5 x 10~ 6 . It begins to dominate over the primary CMB at / ~ 2000. 
The shape of the kinetic SZ power spectrum is similar to the thermal one but the amplitude is 
about 30 times smaller. 
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Fig. 2. — One typical thermal SZ map in our moving mesh hydro (MMH) code simulation. The 
cosmology is a ACDM with r?o = 0.37, Qa = 0.63, $7_b = 0.05, h = 0.7 and as = 1.0. The map size 
is 1.19° x 1.19°. The color represents the SZ temperature fluctuation in the Rayleigh-Jeans regime 
AT/T = —2y. We have omitted the negative sign. The SZ map clearly shows the structures with 
angular scale ~ 1', which is the typical scale of clusters and groups. 



Fig. 3. — The kinetic SZ map in the same simulation at the same field of view as fig. 2. The 
amplitude of temperature fluctuation of the kinetic SZ effect is generally one order of magnitude 
lower than the thermal SZ effect. Gas in the positive AT/T regions is approaching us and those in 
the negative AT/T regions is moving away. 
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Comparing to available CMB observations, our thermal SZ power spectrum is well below the 
upper limit (95% confidence) of ATCA (Subrahmanyan et al. 2000), SUZIE (Church et al. 1997) 
and VLA (Partridge et al. 1997). C\ at I ~ 2000 is consistent with recent indications from the 
CBI experiment (Mason et al. 2001; Sievers et al. 2001) and may suggests a high erg — 1-0. But 
our result is higher than the BIMA l-a result thought it is consistent with the upper limit (95% 
confidence) of the BIMA result(Dawson et al. 2001). We will further discuss these issues below. 

Published theoretical predictions differ a lot in both amplitude and shape (see Springel et al. 
(2001) for a review). Our power spectrum has a higher amplitude than all previous predictions. It is 
also significantly flatter at the range 2000 < / < 15000. The difference in amplitude can be explained 
by the strong dependence of the SZ effect on cosmological parameters, especially as- One expects 
y oc Qsh&s 3 an d Q °c (^_B^) 2fT 8 9 as predicted by various authors (Komatsu & Kitayama 
1999; Seljak, Burwell & Pen 2001; Zhang k Pen 2001). For example, Seljak, Burwell & Pen (2001) 
used the same MMH code with lower resolution (256 3 ) to simulate a universe with as = 0.8. This 
as difference accounts for a factor of 1.6 ~ 2.0 difference in the mean y estimation and 3.8 ~ 7.5 
difference in the C/ estimation. After accounting for these effects, our power spectrum is consistent 
with theirs at small angular scales (I > 3000). But the difference in shape can not be explained in 
this way. For example, even after accounting for the effect of cosmological parameters, our power 
spectrum is still much larger at large angular scales (I ~ 1000) than that of Seljak, Burwell & 
Pen (2001) (at / ~ 1000, ~ 2 times larger). This flatness behavior may be a manifestation of the 
resolution effect. 

We notice that there are numerous high y regions of arcminute or sub-arcminute scales in SZ 
maps as seen in fig. 2. To quantify this phenomenon, we show the dependence of Ng(y > y p ) on 
smoothing angular size (fig. 5). Ng(y > y p ) is a measure of number of structures with angular size 
larger than or comparable to 9. Fig. 5 shows that when the smoothing scale increases from 0.5' to l', 
N(y > y P ) drops significantly in high y regions. This behavior suggests the existence of numerous 
sub-arcminute, high y structures in SZ maps. Higher resolution reveals more such structures. 
Increasing in the number of these objects increases the amplitude of the power spectrum while 
making the power spectrum flatter around the peak. The first effect is obvious and the second 
one can be explained by the Press-Schechter picture. The SZ power spectrum is dominated by 
the halo gas pressure profile fp(r) at all interesting angular scales (Komatsu & Kitayama 1999). 
For a singular isothermal sphere (SIS), f p (r) = n e (r)kBT(r) oc n e (r) oc r~ 2 . Its projection along 
the line of sight then has a 9~ l radial dependence, which produces an angular auto correlation 
function of shape ln(6>). So the resulting power spectrum is nearly flat. A more accurate way to 
see the flatness behavior is to adopt the Limber's equation. The 3D gas pressure power spectrum 
P p (k) oc 5p(k) oc k~ 2 , where S p (k) is the Fourier transform of the gas pressure profile f p (r). From 
Limber's equation, C\ oc J P p (l/x(z),z)f(z)dx °c l~ 2 , thus 1(1 + l)Ci/(2ir) oc 1°. Here, f(z) is the 
redshift dependence of the SZ effect. The halo mass function also plays a role for the flat power 
spectrum. Since SIS profile does not apply to the core of halos, a given cluster is no longer SIS at 
scales smaller than its core size. Smaller clusters take over and extend the flat power spectrum. 




Fig. 4. — The thermal and kinetic SZ effect power spectra in our simulation. They are averaged 
over 40 maps. Dash lines are the corresponding la upper limit and lower limit of the mean 
power spectrum, respectively. Both effects peak at I ~ 10 4 . For the thermal SZ effect, the power 
spectrum is almost flat in the range of 2000 < I < 15000. As a comparison, we show the BIMA 
result (Dawson et al. 2001) under Gaussian assumption. The comparable amplitude between theory 
and observation puts strict constrain on the understanding of the SZ effect. 
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Fig. 5. — N(y > y p ) = J°° p(y)dy is the number of peaks with Compton parameter y bigger than 
y p . The result is averaged over 40 maps. The smoothing function we adopt is the top hat window 
function with radius 8. N(y > y p ) as a function of 6 is a direct measurement of number of structures 
with scale larger than 0. The quick drop of N y>yp from 9 = 0.5' to = l' shows that there are a 
lot of arcminute scale structures. We believe that these structures with large y is responsible for 
the divergence of SZ simulations towards high resolution. 
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The power spectrum in our simulation clearly shows this flatness and suggests the role of these sub- 
arcminute halos. These halos also explain the discrepancy between our simulation and analytical 
predictions, l' corresponds to the comoving size < O.8/1 -1 Mpc for z < 1 where the dominant SZ 
contribution comes from. This physical size corresponds to groups of galaxies. Current predictions 
from the Press-Schechter picture assume a lower mass limit cutoff corresponding to the mass scale 
of groups. In the hierarchical method, the gas window function is a free parameter, which has an 
implicit dependence on this lower mass cutoff (they can be related by the gas density dispersion). 
The cut off for contributions from groups of galaxies in analytical models results in a smaller y and 
power spectrum. 

If only gravitational heating is included, high mass halos and low mass halos should have 
comparable gas fractions. The Press-Schechter formalism predicts many more halos towards the 
low mass end, so one expects pure gravity simulations to have increasing power spectra with 
increasing high resolution. Non-gravitational heating avoids such a divergence. As described in the 
halo model of Pen (1999), non-gravitational heating has two effects. Firstly, the energy injection 
makes the halo gas less clumpy. Then the contribution to the power spectrum at smaller angular 
scales decreases relative to larger angular scales. This could explain the slightly differences between 
our simulation, CBI and BIMA results. Secondly, the non-gravitational energy injection increases 
the thermal energy of the gas. For halos with mass lower than some threshold, the gravity can not 
hold gas and most of the gas is ejected from these halos, as must have been the case for galactic size 
halos. This provide a reasonable lower mass limit cut off in the Press-Schechter picture. Thus the 
amplitude and shape of the SZ power spectrum is a sensitive measurement of the non-gravitational 
energy injection. To obtain a better understanding of the SZ effect, other physicical processes such 
as radiative cooling need to be considered. Since radiative cooling through thermal bremsstrahlung 
is a p 2 process, it becomes relevant only at scales < 100 kpc, which corresponds to / S> 10 . 
These angular scales are not observable by any planed experiments, so we neglect the discussion of 
radiative cooling in this paper. Our current simulation has reached the resolution needed to see the 
contribution from small halos and the predicted SZ amplitude is already near the observed values. 
We expected to be able to observe the effects of non-gravitational heating from galaxy winds and 
other sources in the upcoming experiments. 

In order to solve the discrepancy problem completely, differences between codes must be con- 
sidered. Our group is currently running different codes with identical initial condition, identical 
cosmological parameters and various resolutions from 64 3 to 512 3 (Codes Comparison Program 
2002). 

Our prediction is consistent with the BIMA 95% confidence result. But the comparison between 
observations and theoretical predictions needs further investigation. On the simulation part, as we 
discuss above, the detailed normalization depends on <7g, non-gravitational heating and resolution 
effects. Observations will put strong constrains on these aspects. On the observational part, the 
BIMA result needs to be reconsidered. It covers the sky where no strong SZ temperature distortion 
due to known galaxy clusters exists. Since clusters contribute a significant if not dominant fraction 
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to the SZ power spectrum, the BIMA result may be smaller than the statistical mean value. 
Furthermore, the error of the BIMA result is estimated under the Gaussian assumption, but as we 
will see below, the SZ effect is highly non-Gaussian around the BIMA central multipole I = 5530. 
The strong non-Gaussianity increases the intrinsic error of the SZ power spectrum measurement. 
According to fig. 6, the error in the power spectrum measurement caused by the SZ effect is about 

3 times larger than the corresponding Gaussian case. 

Nontheless, these results demonstrate the convergence between theory and observation. In the 
near future, routine and accurate measurement of the SZ effect will be possible in random fields. 
It will put stronger requirement for our theoretical understanding of the SZ effect and allows us 
to study the effects of non-gravitational heating, radiative cooling and the thermal history of the 
IGM. 

3.2. The SZ non-Gaussianity 

In contrast to the primary CMB, the SZ effect is non-Gaussian, arising from the non-linearity 
of the intervening gas. This non-Gaussianity affects the error analysis and may help to separate 
the SZ effect from the primary CMB in observation. To quantify these effects, we smooth SZ maps 
using a top hat window of radius 9 and measure the kurtosis of the smoothed y. The kurtosis 

4 = ( y ~p 3 is generally a function of smoothing scale. For a Gaussian signal, ©4 = 0. Fig. 

a y 

6 shows that O4 ~ 200 at small scales as 9 — > 0. This result is consistent with the prediction 
from our hierarchical model approach of the SZ effect (Zhang & Pen 2001). We predicted, at small 
angular scales, 64 ~ S\(T^{z ~ 1) ~ 200. Here, 64 = -^p- ~ 40 is a hierarchical model coefficient 
(Scoccimarro & Frieman 1999) and a g (z ~ 1) ~ 5 is the gas density dispersion at z ~ 1. At large 
angular scales 9 ~ 20', 84 3> 1 and reflects the strong non-Gaussianity at this scale. It means 
that, at this angular scale, the dominant contribution to the y parameter is from highly non-linear 
regions and there is still strong correlation at angular scales down to I ~ 1000, as can be seen from 
the SZ power spectrum. 

The kurtosis in multipole space 04 = (|a| 4 )/(|o| 2 ) 2 — 3 . As usual, a = a\ m defined by 
0(n) = Ylim a im^i m (n) are multipole modes of and C\ = Y^-i \ a im\ 2 /{^ + 1) = (|«| 2 )- Here, 
yj m are the spherical harmonics. (...) is averaged over all m = —I,..., I and all maps. We 
calculate indirectly through the map-map variance 0m(O of Ci, which is related to 04 by 
cm(0 = Q\/(a4 + 2)/[(2Z + l)AZ/ map ]. / map is the fractional sky coverage of each map which 
reflects the cosmic variance and Al = 27r/# map is the I bin size in our grid map with map = 1.19°. 
The factor 21 + 1 arises because C\ is averaged over 2^ + 1 independent a\ m modes. In our calcula- 
tion, C\ of each map is obtained using FFT under the flat sky approximation. Then we obtain the 
map-map variance cjm(0- O ne might expect 04 to have similar behavior as 04 at corresponding 
scales, for example, at small angular scales (large I), 04 » 2. But fig. 6 shows the opposite. 
04 ~ 50 around I ~ 1000, rises slowly until I ~ 3000 and approaches Gaussian (04 = 0) quickly 
after that. In fact, ©4 and 04 cannot be directly compared. The multipole modes are non-local, 
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Fig. 6. — The non-Gaussianity of the SZ effects in real space (left panel) and multipole space (right 
panel). In real space, we measure the skewness 03 and kurtosis ©4. For symmetrically distributed 
signals, O3 = 0, as in the kinetic SZ case. For a Gaussian signal, ©4 = 0. The left panel shows that 
the thermal SZ effect is highly non-Gaussian even to the angular size 6 ~ 0.3°. The kinetic SZ is 
mildly non-Gaussian and approaches Gaussian soon towards large angular scales. The right panel 
shows the non-Gaussianity in multipole space. 04 is the kurtosis in the multipole space. Because 
multiple modes are not local and averaged over many uncorrelated patches, the corresponding 
non-Gaussianity is much smaller than in the real space. 
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Fig. 7. — The log 10 y PDF P(log 10 y). We use the top hat window with radius 9 to smooth the SZ 
map. P(log 10 y) is roughly Gaussian. So the y PDF P(y) = P(log 10 y)/(y In 10) is non-Gaussian 
even with a large smoothing scale as shown in fig. 6. 
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and averaged over many patches. If these patches are spatially uncorrelated, then the central limit 
theorem assures the multipole modes to be Gaussian independent of the actual non-Gaussianity 
of the patches. From this viewpoint, we can estimate the relation between 04 and ©4. Yj ra (h) 
are quickly fluctuating functions with period A6 ~ 2-ir/l. Two signals separated by a distance 
larger than 4A6> will have effectively no correlation. So, for a I mode, a; m is approximately the 
sum of N ~ #^ ap /(4A#) 2 uncorrelated patches. The non-Gaussianity of each patch is roughly 
characterized by 64. Then, 04 ~ ©4/iV, which may explain the behavior of 04. 

The shape of a±(l) is similar to the shape of the pressure bias b p (k,z) (fig. 2, Zhang & Pen 
(2001)) defined by bp(k, z) = Pp(k, z)/Ps(k, z) where P p and P$ are the corresponding gas pressure 
power spectrum and dark matter density power spectrum. This similarity suggests a common 
origin. In the hierarchical approach, a±{l) is related to b p (k, z) and we expect 0,4(1) to keep the 
similar shape but with smaller amplitude variation than b p , as illustrated by the behavior of a^l) 
in our simulation and b p (k, z) in our analytical model. 

We further show the y PDF (fig. 7) as a function of smoothing scale. The distribution of 

lo 

gio V i s roughly Gaussian, especially for large smoothing scales, as suggested by Seljak, Burwell 
k. Pen (2001). This behavior reflects the non-Gaussianity of the y parameter, up to angular scales 
~ 10'. The distribution of y is asymmetric since the distribution of log 10 y is nearly symmetric. We 
quantify it by the skewness ©3 = ((y — y) 3 )/o"^ as a function of smoothing scale. For symmetric 
y distribution, ©3 = 0. Fig. 6 shows that ©3 ~ 10 when 9^0. This result again agrees with 
our hierarchical method prediction. We predicted, at small angular scales, ©3 ~ Ss(Jg' 5 (z ~ 1) ~ 
10 where S3 = jjph ~ 5 (Scoccimarro & Frieman 1999). The positive sign of ©3 reflects the 
numerousness of high y regions, as suggested in fig. 5. 

The kinetic SZ effect is determined by the gas momentum along lines of sight. The fractional 
temperature change @(n) = —ot Jq C n e (ln)™dl. Since the direction of the gas velocity v is 
random, we would expect ©3 = for a sufficiently large sky. The contribution from non-linear 
structures is partly cancelled out and thus we expect mildly non-Gaussianity. As shown in fig. (6), 
the kinetic SZ ©4 and 04 are significantly smaller than the corresponding thermal SZ ©4 and 04 
(fig. 6). The kinetic SZ effect is nearly Gaussian at angular scales ~ 20'. 

4. Simulated observations of the SZ effect 

In real SZ observations, instrumental noise and primary CMB cause additional errors in the SZ 
statistics such as the power spectrum and peak number counts. We need to estimate these effects 
to derive optimal observing strategies to measure these statistics in the presence of noise. With 
these observational errors, our methods (§2) to extract 3D gas information is limited and we must 
check their feasibility. In this section, we take AMIBA as our target to address these issues. 
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4.1. AMIBA 

AMIBA is a 19 element interferometer with 1.2 meter dishes. All dishes are closely packed in 
three concentric rings. It operates at f CC ntcr = 90 Ghz with Au = 16 Ghz, system noise T sys = 100 K 
and system efficiency rf ~ 0.7. At this frequency, ~ — 1.6y with i>T(90Ghz) ~ 0.8. The goal of this 
experiment is to image maps of the CMB sky with arc minute resolution. We consider observations 
with fixed integration time and aim at finding the optimal sky area O and sky fractional coverage 
/sky = ^/47r for a given statistics. 

For closely packed interferometers observing such weak signals, the ground fringe can be a 
major source of noise. To eliminate the ground fringe, AMIBA plans to drift scan: the telescope 
is parked while the sky drifts by. This assures that the ground fringe remains constant with time. 
The mean value of each fringe is then subtracted from the scan, cleanly eliminating the ground. 
A field is mosaicked by a series of adjacent scans, and the most uniform coverage is achieved by 
incrementally offsetting the pointing center on each scan to yield a finely sampled 2-D map. The 
raw output of the experiment are correlations, two for each baseline, polarization and frequency 
channel, corresponding to the real and imaginary correlator outputs. We can think of each of 
these outputs to correspond to an image of the sky filtered through some anisotropic beam. As 
a first step, we can combine degenerate baselines and polarizations, reducing the 171 baselines to 
30 non-degenerate baselines. Since the CMB is expected to not be significantly unpolarized, we 
can combine the two polarization channels, leading to 60 raw maps per frequency channel. These 
maps can be merged optimally into one global map by convolving each map with its own beam, 
and scaling each map to the same noise level, and coadding these maps, resulting in a 'clean' map. 
Each of the constituent maps had explicitly white noise, so the noise statistics of the summed map 
are easily computable. The 'clean' map can be deconvolved by the natural beam, resulting in a 
'natural' map. This 'natural' map is an image of the sky convolved with the natural beam of the 
telescope plus white noise. The angle averaged natural beam is shown in fig. 9. The CMB intensity 
fluctuation 8I V measured by AMIBA has three components: the primary CMB, the SZ effect and 
the instrumental noise. It is related to the temperature fluctuation by 



5I V xexp(x) 



5T\ /5T\ 5T N 



/CMB exp(x) - 1 l\T J CMB V T J sz T RJ 



(6) 



Here, (^) CMB and (§-) sz are the corresponding temperature fluctuations seen through the AMIBA 
natural beam, x = hu/ksTcMB- We have chosen the normalization of the beam such that the noise 
power spectrum is white with 

C N = 47rT s 2 ys / sky /(2A^7 ? 2 )/T2 J . (7) 

Trj = Xcmb^ 2 exp(x)/[exp(x) — 1] = 2.22 K is the Raleigh-Jeans equivalent CMB temperature 
at v = 90 Ghz. The factor of 2 is due to the two polarizations of the AMIBA experiment. For a 
single dish of infinite aperture, with a single pixel detector, the window would be identically one: 
the scanned image is just the CMB distribution on the sky. Since AMIBA has many detectors, 
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one can combine them to either lower the noise, or to boost the signal. We have chosen to use 
the equivalent noise of a single pixel single dish experiment, and normalized the beam accordingly. 
We calculate the natural beam Wn(1) by equation (17) and (46) of Pen et al. (2002). The natural 
beam in multipole space and real space are shown in fig. 9. The natural beam has a FWHM of 
2'. It peaks at Zj ~ 27rA/Z)j. Here, A ~ 3.3 mm is the AMIBA operating wavelength and Di is 
the distance of the ith baselines. The first peak l\ eak ~ 2273 corresponds to the shortest base line 
D\ = 1.2 m. At this angular scale, the sum of all the baselines improves throughput by a factor of 
almost three (e ~ 2.77). This is analogous to having a nine pixel detector on a single dish. More 
detailed definitions of the beams and strategies are described in detail in (Pen et al. 2002). 

In our simulated pipeline, we add primary CMB map fluctuations to our simulated sky maps 
using CMBFAST-generated (Seljak & Zaldarriaga 1996) power spectra with same cosmological 
parameters except for the use of COBE normalization for eg, which is slightly different from our 
simulation value. As we will show below, the primary CMB is not the main source of noise for the 
SZ power spectrum measurement at AMIBA angular scales (/ > 2000) and is negligible in SZ cluster 
searches. So the effect of this inconsistent as is insignificant in our analysis. Adding SZ maps and 
CMB maps, we obtain simulated sky maps (fig. 8). In these maps, SZ structures, especially those 
caused by diffuse IGM are superimposed with the primary CMB. We then convolve these maps with 
the natural beam. The beam function decreases quickly to zero towards large angular scales where 
the primary CMB dominates, so it efficiently filters most primary CMB structures larger than this 
scale. We then add the noise given by eqn. (7) to this map. The normalization in the beaming and 
filtering (described below) process is arbitrary. We choose the normalization in such a way that 
the power spectrum of the map at the scale of the peak response does not change after beaming 
or filtering. It corresponds to normalizing the global maximum of the beam function to be unity. 
Under such normalization, when we add instrumental noise to the simulated CMB+SZ map, the 
noise is depressed by a normalization factor e ~ 2.77. Our 1.19° map has 2048 2 pixels, so the white 
instrumental noise dominates on small scales. For a 20 hours/deg 2 survey, the dispersion of the noise 
temperature fluctuation <tat = T sys / yj (2Az/i pixc i/ (r/Tjjje) ~ 0.002 3> y. t pixc \ is the observing time 
for each pixel. So all signals are hidden under the instrumental noise. One needs further filtering to 
obtain an image that is not overwhelmed by the small scale noise. Ignoring the CMB fluctuations, 
we know that point sources have the shape of the beam. A point sources optimized search would 
convolve the natural map with the shape of the beam, and peaks in this map correspond to the 
maximum likelihood locations of point sources. We can think of the CMB as a further source of 
noise, and filter that away as well. In this case, the filter depends on the ratio of CMB and noise 
amplitudes, i.e. the actual integration time. One can also optimize a filter for a structure of a 
known intrinsic shape. Since clusters are approximately isothermal, a better filter would be one 
matched to the extended isothermal nature, where objects in the natural map have the shape of an 
isothermal sphere convolved with the natural beam. Combining all these consideration, the optimal 
filter for noise cleaning is given by W c (l) = W F (l)W N (l)/ '(1 + W£C, CMB /C N ) (fig. 9). W F {1) is the 
Fourier transform of the source intrinsic shape Wf(r). For point sources, Wf(1) = 1. For clusters, 
on scales smaller than the cluster virial radius and larger than the core radius, y(6) oc _1 (9 is 




Fig. 8. — A combined SZ+CMB map. At large angular scales, the primary CMB dominates and 
smears the structure of the thermal SZ effect. 
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the angular distance to the cluster center). Then Wf{1) oc l^ 1 . Since the corresponding angular 
size of the core radius is much smaller than the beam size, the above approximation is sufficiently 
accurate. We shown the filter with 280 hours per square degree scan rate in figure 9. It has a 
FWHM of 2 . In multipole space, it drops to near zero at angular scales I ~ 1500 and I > 9000 
and peak at cluster scales I ~ 3000. So it is efficient to filter away the residual primary CMB and 
the instrumental noise while amplifying SZ signals. We tried different sky scan rates to find the 
optimal survey strategy. We show final resulting maps in fig. (10) with 280 hours per square degree 
scan rate. With the simulations we have the luxury of seeing maps with (right panels of fig. 10) 
and without noise (left panels of fig. 10). We discuss the simulated AMIBA measurement of the 
SZ power spectrum and cluster searching in the next two subsections. 



4.2. The SZ power spectrum in the simulated AMIBA experiment 

One of the key goals of AMIBA is to measure the SZ power spectrum. As discussed in §3, 
it is a sensitive measure of the gas thermal history. Furthermore, combining the cross-correlation 
with photometric galaxy surveys, an SZ survey can measure the IGM pressure power spectrum as 
a function of redshift. This gives us access to the evolution of the IGM state. For the purpose 
of the SZ power spectrum estimation, the two noise sources, the primary CMB and the thermal 
instrument noise, are both Gaussian. The intrinsic SZ variance causes further error. The combined 
error for the power spectrum estimation is 



-ySZ 



Acf z = I °4(QCg z (Q + 2 [Csz(Q + Cp MB + CljWh _ 



(21 + 1)AZ/ Sky 

If the SZ effect is Gaussian (04 = 0), we recover the usual expression of the error. We take a large 
bin width Al = 1/4 to estimate the error. Since the SZ power spectrum is nearly flat in the range 
I ~ 2000 — 15000, this choice of bin size does not lose significant information. At large angular scales 
(I < 1500), the error from primary CMB dominates, and at small scales (I > 5000), the instrumental 
noise dominates. In the intermediate scales, the intrinsic error of the SZ effect dominates. Different 
errors have different dependences on the sky coverage. Increasing /f s k y while fixing the integration 
time, the errors from primary CMB and SZ effect decrease but the the one from the instrumental 
noise increases, so there exists a optimal sky coverage for a given integration time. Since the errors 
from primary CMB and SZ effect are both scale as f s ^ 2 and the instrumental noise scales as f^[y, 
the minimum error is obtained when (a 4 + 2)Cf z + 2C, CMB + 4C S zC CMB = 2C%/W^, which gives 
the optimal sky coverage. For 1000 hours of observing, several hundred square degree sky coverage 
is nearly optimal (fig. 12). Fig. 11 shows that AMIBA is able to measure the SZ power spectrum 
with an accuracy of ~ 40% for I between 2000 and 5000 in 1000 hour observing of 100 square 
degrees of sky. The strong dependence of the SZ effect on as strongly affects our error estimation. 
A smaller as reduces the signal and S/N. In order to keep S/N, smaller sky coverage is needed to 
reduce the noise. The dependence of the optimal sky coverage f°£ y on as is shown in fig. 12. In 
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Fig. 9. — The natural beam and optimal filter functions of AMIBA in multipole (left panel) and 
real space (right panel). The main goal of the beaming and filtering is to filter away the primary 
CMB, which dominates at large scale (I < 1000) and instrumental noise at small scale(/ — ► oo). 
This goal is clearly illustrated in the large and small / behavior of these functions. The optimal 
filter depends on the noise amplitude. We show the case for AMIBA scan rate 280 hours deg -2 . 
This filter has a ~ i FWHM. 




Fig. 10. — Filtered maps without noise (left panel) and with noise (right panel) of AMIBA 390 
hours observing. Many peaks in the total map are not real signal peaks while some peaks in the 
clean map disappear in the total map. For AMIBA's frequency, = — 1.6y. For clarity, we plot 
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Fig. 11. — Errors in the AMIBA measurement of SZ power spectrum (left panel) and cross cor- 
relation (right panel) with SDSS galaxies. The SZ power spectrum is 0.8 2 times of the one in fig. 
4 since at AMIBA operating frequency v = 90 Ghz, = — 1.6y as contrast to = — 2y in the 
Rayleigh-Jeans regime . The thick solid error bar is the Gaussian variance, the thin solid error bar 
is the actual variance calculated from a 40 map ensemble and the dot error bar is the total error 
including primary CMB and instrumental noise. For a 1000 hours survey of a 100 deg 2 area, the 
accuracy in the SZ power spectrum measurement is about 40% at the range of 2000 < I < 5000. 
In the same survey, the measured cross correlation has about 20% accuracy at similar I range. 
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Fig. 12. — The dependence of the optimal sky coverage in the SZ power spectrum measurement on 
<7g in a 1000 hour AMIBA survey. Our simulation used as = 1.0. The as = 0.9 and a% = 0.8 cases 
are estimated assuming the strongest dependence of C[ on a$: C\ oc a\. 
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this estimation, we have assumed the most extreme dependence of the SZ effect on as, namely, 
C S z{l) OC of. 

If we cross correlate the observed SZ effect with SDSS, we can extract the underlying 3D gas 
pressure power spectrum and pressure-galaxy cross correlation. To test its feasibility, we estimate 
the error in the angular cross correlation measurement by 

AC SZ - G (Z) _ l l + r~2(l + C^/Csz + C s yCsz)(l + C^/QG) 
C^ G (l) V (2Z + l)A// sky • W 

SDSS will cover / sk y = 1/4 of the sky and will detect Nq — 5 x 10 7 galaxies with photometry in 
five bands. The Poisson noise in SDSS has the power spectrum C% = Airf^ ss /N G . We assume a 
linear bias between galaxy number overdensity and dark matter overdensity to calculate the galaxy 
surface density power spectrum C G (l). Since SDSS is flux limited, we take the galaxy selection 
function %± = 3z 2 /2/(z m /lA12) 3 exp[-(lA12z/z m ) 3 / 2 ] (Baugh & Efstathiou 1993) with a SDSS 
fit 

z m — 0.33 (Dodelson et al. 2001). The galaxy-SZ power spectrum C^^' G (l) is estimated by 
the SZ-galaxy cross correlation coefficient r = C SZ ' G / V C SZ C G . We choose r = 0.7 (Zhang & Pen 
2001) as predicted by the hierarchical model. We show the error in fig. 11. The larger the sky 
coverage, the smaller the error. With the optimal scan rate for SZ power spectrum measurement, 
the accuracy in the cross correlation measurement is about 20%. 



4.3. AMIBA cluster search 



Another important goal of AMIBA is searching for clusters. In the power spectrum mea- 
surement, the observable is the direct sum of signal and noise, so the noise contribution can be 
subtracted linearly in the power spectrum. The only net effect is an increase in the error bars, which 
is important only at very small angular scales or very large scales. But when counting peaks in a 
map, effects of noise are much more complicated and are not easily interpreted. Noise introduces 
false peaks in the observed SZ maps, changes the value of real peaks, shifts the peak positions and 
even makes real peaks disappear. So noise affects the measurement of cluster counts, richness and 
positions. Thus we would want a much longer integration time for the purpose of cluster search. 
We count peaks in the filtered maps with and without noise and quantify these effects as follows. 
(1) We estimate the accuracy of y measurements by calculating the y dispersion of noise . (2) 
We distinguish real peaks from false peaks by comparing the position and amplitude of each peak 
in the clean SZ maps without noises and in the total maps with noises (fig. 13). One FWHM is 
roughly the size of noise structures after filtering and corresponds to the maximum position shift 
noise can exert to real peaks. 2(jy is roughly the maximum peak amplitude change that noise can 
cause. So, if a peak in a total map whose distance to the nearest peak in the corresponding clean 
SZ map is less than one FWHM of the filter and its value is in the 20y range of the real value, 
we classify it as real. (3) We estimate the accuracy of y peak CDF by comparing clean SZ maps, 
noise maps and total maps (fig. 14). For signal peaks with y » <Jy , the signal peaks remain to 
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be peaks in the total map. Noise mainly changes the amplitude of the zth real signal from yi to 
Up = Ui + y N i n the total map. Here, y N is the value of the noise and y p is the value in the total 
map. Since noise is Gaussian, we know the distribution of y N , which can be described by P(y, &y), 
the probability for the Gaussian noise with dispersion <7y to have value bigger than y. Then we 
can related the CDF of peaks in clean maps to the one of total maps by 

N 

N tot (y >y P ) = Yl p (y P - ^ O- ( 10 ) 

i 

This relation gives a good fit in the y » regime. For signals with y ~ a^, noise changes both y 
value and positions of some peaks, makes some peaks disappear and introduces a large fraction of 
false peaks. The sum of N tot {y > y p ) in eqn. (10) and N no i se (y > y p ) considers the effect that noise 
introduces false peaks and we expect that it would give a good fit in the region where y ~ 2<jy . We 
show the modelled CDF of peaks in total maps following above procedures in fig. 15. The result is 
well fitted in the range y > 2<jy with better than 30% accuracy. 

The optimal survey should find clusters as quickly as possible and measure the y parameter 
and peak CDF as accurately as possible. We find that the optimal scan rate for cluster searching is 
about 280 hours deg -2 . At this rate, the rate of cluster detection is about 1 every 7 hours allowing 
for a false positive rate 20%. The measured cluster CDF N(y > y p ) is accurate to 30% up to 
N ~ 60 deg -2 . We recall that the optimal scan rate for SZ power spectrum measurement is about 
10 hours deg -2 . At this rate, the cluster detection rate is about 1 cluster every 60 hours. We can 
consider the effect of cosmology on our estimation. A smaller a$ reduces signals and therefore the 
cluster detection rate. We assume the extremest dependence of the y parameter on as, namely, 
y oc erf as predicted by Zhang & Pen (2001). For ag = 0.9, the detection rate decreases by a factor 
of four. But the optimal scan rate remains approximately the same. 

At this optimal scan rate, in a 1000 hour AMIBA survey, several hundred clusters can be found. 
Comparing the SZ cluster counts to Press-Schechter predictions or X-ray surveys (§2) allows us to 
reconstruct the thermal history of the IGM and non-gravitational heating, which is expected to 
arise from galaxy formation feedback (Pen 1999). 

5. Conclusion 

We have performed the largest high resolution SZ simulations to date, and analysed the results. 
We found further increases in the power spectrum relative to previous simulations, and significantly 
more small scale structures. This trend is confirmed in Press-Schechter estimates, and suggests 
nominally increasing power spectra. The small structure behavior may not be a robust prediction, 
since non-gravitational effects will significantly modify those scales. We have examined the skewness 
and kurtosis on the sky maps, found a strong non-Gaussianity on subdegree scales and confirmed 
the log- normal distribution found in previous studies. The Gaussian estimates of power spectrum 
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sample variances are less severely affected, due to the averaging of many patches by each Fourier 
modes. But its effect on the power spectrum error analysis is significant. 

We simulated SZ observations with AMIBA, and analyzed the sensitivity for different scan 
rates. In a 1000 hour survey, the optimal strategy for power spectrum estimations is to scan 
several hundred square degrees. The SZ angular power spectrum measured in such survey can be 
determined to an accuracy of ~ 40% over a range of 2000 < I < 5000. A cross correlation with 
SDSS should allow an accuracy of 20% in the cross correlation measurement, which suggests that 
the time resolved measurement of the pressure power spectrum is then possible. This scan rate 
results in a low detection rate of clusters of galaxies, approximately one per 60 hours with a false 
positive rate of 20%. For the purpose of cluster search, the optimal scan rate is around 280 hours 
per square degree, which could find 1 cluster every 7 hours, with the accuracy of 30% in the cluster 
y measurement and 30% in the cluster CDF up to N(y > y p ) ~ 60 deg~ 2 . 

The predicted SZ power spectrum is consistent with recent indications from the CBI experiment 
(Mason et al. 2001; Sievers et al. 2001) and the BIMA upper limit (95% confidence). But it is higher 
than the BIMA 1-cr result. This may be a first indication of IGM non-gravitational feedback. Future 
blank sky surveys with data analysis considering actual SZ non-Gaussianity will provide us with a 
quantitative understanding of the thermal history of the universe. 
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Fig. 13. — The fraction of real peaks as a function of number of peaks in simulated AMIBA 
observing (left panel) with a% = 1.0 and the expected result scaled for a different a% = 0.9 (right 
panel). In our simulation, the optimal cluster searching rate is ~ 1 cluster every 7 hours, which 
drops to ~ 1 cluster every 30 hours for erg = 0.9. We have assumed the extremest y dependence on 
<7 8 : y oc erf. 
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Fig. 14. — Peak distribution N{y > y p ) in simulated AMIBA observing. All results are averaged 
from 40 maps, the 'TSZ' line is the pure SZ peak distribution without noise and the 'total' line is 
what would be measured in our simulated AMIBA experiment. We see that when y < 2<r^ , peaks 
due to noise dominate the peak CDF and makes the measurement of peak CDF unreliable at this 
region. 
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Fig. 15. — Modelled N(y > y p ). The line 'model 1' is calculated from eqn. (10), which consider 
the effect that noise changes the amplitude of peaks. The line 'model 2' is the sum of 'model 1' 
and noise peak CDF, which considers the effect that noise introduces false peaks. This result fits 
the CDF of peaks in the total maps down to y p ~ 2<jy with better than 30% accuracy. 



